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Abstract 

In a multi-type totally asymmetric simple exclusion process (TASEP) on the line, each site of Z is 
occupied by a particle labeled with a number and two neighboring particles are interchanged at rate 
one if their labels are in increasing order. Consider the process with the initial configuration where 
each particle is labeled by its position. It is known that in this case a.s. each particle has an asymptotic 
speed which is distributed uniformly on [-1, 1]. We study the joint distribution of these speeds: the 
TASEP speed process. 

We prove that the TASEP speed process is stationary with respect to the multi-type TASEP dy- 
namics. Consequently, every ergodic stationary measure is given as a projection of the speed process 
measure. 

By relating this form to the known stationary measures for multi-type TASEPs with finitely many 
types we compute several marginals of the speed process, including the joint density of two and three 
consecutive speeds. One striking property of the distribution is that two speeds are equal with positive 
probability and for any given particle there are infinitely many others with the same speed. 

We also study the (partially) asymmetric simple exclusion process (ASEP). We prove that the 
ASEP with the above initial configuration has a certain symmetry. This allows us to extend some of 
our results, including the stationarity and description of all ergodic stationary measures, also to the 
ASEP 

1 Introduction 

The exclusion process on a graph describes a system of particles performing continuous time random 
walks, interacting with other particles via exclusion: attempted jumps to occupied sites are suppressed. 
When the graph is Z and particles jump only to the right at rate one the process is called the totally 
asymmetric simple exclusion process (TASEP). We denote configurations with rj e {I, oo}^ where particles 
are denoted by 1 and empty sites by ooj^ The TASEP is a Markov process with generator 

Lm = Y,f((^nri)-m (1) 

n 
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where cr„ is the operation that sorts the coordinates at «, « + 1 in decreasing order: 



(o-„77)„ = 



max(?7„,?7„+i), 



(cr„77)„+i = min(77,„?7„+i), 



(cr„77)i = r]k if kin, n + I. (2) 



A second class particle is an extra particle in the system trying to perform the same random walk 
while being treated by the normal (first class) particles as an empty site. It is an intermediate state 
between a particle and an empty site, and is denoted by a 2|1 This means that the second class particle 
will jump to the left if there is a first class particle there who decides to jump onto the second class 
particle. This is still a Markov process, with the same generator ^ and state space {1,2, oo}^. Note that 
empty sites can just be considered as particles with the highest possible class. Thus we can equally well 
consider state space {1, 2, 3}^ with holes represented by 3s. 

More generally, we shall consider the multi-type TASEP which has the same generator with state 
space R^. Thus we allow particle classes to be non-integers or negative numbers. If there are particles 
with maximal class they can be considered to be holes. A special case is the A^-type TASEP (without 
holes) where all particles have classes in { 1 , . . . , A'^} . If particles of class N are interpreted as holes instead 
of maximally classed particles, this process becomes the traditional (A^ - l)-type TASEP (with holes). To 
avoid confusion, from here on all multi-type configurations shall be without holes. (Holes will appear 
only in individual lines in the multi-line configurations defined below.) 

The following result is this paper's foundation. We let Y{t) denote the TASEP configuration at time 
t, with Y„{t) the value at position n. This strengthens results of Ferrari and Kipnis [7] that get the same 
limit in distribution. 

Theorem 1.1 (Mountford-Guiol, [13]). Consider the TASEP with initial condition 



where U is a uniform random variable on[-\,\]. 

Thus a second class particle with first class particles to its left and third class particles to its right 
"chooses" a speed U, uniform in [-1, 1] and follows that speed: X(t) ~ Ut. (See [9, 10] for alternative 
proofs of [Theorem 1.11 ) 

Now, consider any other starting configuration such that y„(0) < ^0(0) for all n < and Y„{0) > Yq{0) 
for all n > 0. The particle starting at does not distinguish between higher classes, or between lower 
classes, so its trajectory has the same law. This applies in particular to every particle in a multi-type 
TASEP Y with starting configuration y„(0) = n. Let Xn{t) be the location of particle n at time t, so that 
Yx„{t){t) = n (^(0 is the inverse permutation of Y{t)). An immediate consequence is the following: 



1 n < 0, 
y„(0)= 2 n = Q, 
3 n > 0. 



LetXit) denote the position of the second class particle at time t, defined by Yx{t){t) = 2. Then 



a.s. 



^ u, 



^kxh class particles will be denoted by k, even for k -0. That is why it is convenient to use oo for holes rather than 0. 
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Figure 1: The speed process: Simulation of Un for I < n < 5000, from a simulation run to time 700000. 

Corollary 1.2 (The speed process). In the TASEP with starting configuration Yn(0) = n, a.s. every parti- 
cle has a speed: For every n 

X„{t)-n a.s. 

> Un, 

t t^oo 

where {Un}nez is a family of random variables, each uniform on [-1, 1]. 

Definition 1.3. The process {Un}nez is called the TASEP speed process. Its distribution is denoted by fi. 

Thus // is a measure supported on [-1, 1]^. It is clear from simulations (and our results below) that fi 
is not a product measure, i.e. that the speeds are not independent. Figure [U shows a portion of the process. 
Some aspects of this process were studied in [6]. 

1.1 Main results 

In order to study the TASEP speed process we prove two results, which are our main tools in understand- 
ing the joint distributions of speeds. These results are of significant interest in and of themselves. The 
following is a new and surprising symmetry of the TASEP. A version of this theorem was proved in [2], in 
the context of the TASEP on finite intervals. We extend it here also to the asefI (defined in Section [T31) . 

Theorem 1.4. Consider the starting configuration Y„(0) = n and X„(t) as above. For any fixed t > the 
process {X„(t)]nez has the same distribution as {Y„{t)}nez- This holds also for the ASEP. 

At any time t we have that X{t) and Y{t) are permutations of Z, one the inverse of the other. Thus this 
theorem implies that Y{t) as a permutation has the same law as its inverse. It is not hard to see that this 
holds only for a fixed t, and not as processes in t (e.g. Xo(0 changes by at most 1 at each jump). 

^Some sources use PASEP/ASEP respectively for what other sources call ASEP/TASEP (PASEP stands for Partially. . . ). 
We adopt the latter convention. 
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Figure 2: The joint distribution of Uq, U i : Based on 5000 pairs from a simulation run to time 25000. 

The next resuh gives additional motivation for considering the speed process, as it relates its law jd to 
stationary measures of the multi-type TASEP (and ASEP). 

Theorem 1.5. fi is itself a stationary measure for the TASEP: the unique ergodic stationary measure 
which has marginals uniform on [-1,1]. 

This means that if we consider a TASEP with initial configuration 7(0) in [-1, 1]^ with law fi then at 
any time t the distribution of Y{t) is also given by /i. 

It is known that the N-type process has ergodic stationary measures, and that the distribution of ¥„ 
among the classes determines this distribution uniquely. Standard techniques (see below) can be used to 
show that the same holds also with infinitely many classes. Specifically, for any distribution on R there is 
a unique ergodic stationary measure for the TASEP with Yq (and any Y„) having that distribution. For any 
two non-atomic distributions on R, these measures are related by applying pointwise a non-decreasing 
function to the particle classes (see lLemma 5.31) . so every such measure can be deduced from the measure 
with marginals uniform on [-1, 1]. If a distribution has atoms then the corresponding stationary measure 
can still be deduced from the speed process' law fi in the same way, but the operation is non-reversible. 
Thus we have the following characterization: 

Corollary 1.6. Every ergodic stationary measure for the TASEP can be deduced from n by taking the 
law of{F(Un)}nezfor some non-decreasing function F : [-1, 1] ^ R. 

1.2 Results: joint distribution 

Computer simulations suggested early on that Uq, U\ are not independent (see Figure |2l). Recent results 
of Ferrari, Goncalves and Martin [6] confirm this prediction. They proved (among other things) that the 
probability that particle eventually overtakes particle 1 (we identify a particle with its class) is 2/3. It 
follows that P(i7o ^ Ui) > 2/3 (not necessarily equal since Uq = Ui does not a priori imply overtaking). 
Our first theorem describing the joint distribution of speeds is the following: 



The TASEP speed process 



Amir, Angel & Valko 



5 



Theorem 1.7. The joint distribution of(Uo, U[), supported on [-1, 1]^, is 



f{x,y) dxdy + g{x)l{x = y)dx 



with 



I 



4 ^>y^ 
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In particular, V(Uo > C/i) = 1/2, V(Uo = C/i) = 1/6 andF(Uo < t/i) = 1/3. 

Remarks Note that the density in {Uq < Ui} (linear in [/i - Uq, so that there is repulsion between 
the speeds) can be deduced using only ITheorem 1 .41 (we do not include this argument here). However, 
proving the — seemingly simpler — constant density on{Uo > Ui) and deriving the singular component 
on the diagonal requires the power of ITheorem 1.51 It is interesting to compare the power of ITheorem 1.41 
with that of the methods of [6] . It appears that both methods run into similar difficulties and have similar 
consequences, suggesting a fundamental connection (there are also some parallels in the proofs). Specif- 
ically, can the density in the region {Uq < Ui) can be derived using the techniques of [9]? Finally, it is 
interesting that our proof relies non-trivially on the extension of the TASEP to infinitely many different 
classes of particles, though the question and answer can both be posed using only 4 classes (including 
holes). A similar remark holds about some other results below as well. 

Additional information about the joint distribution of speeds is derived in Section|71 We derive certain 
properties of the n-dimensional marginals of p and in lTheorem 7.71 we compute the joint distribution of 
three consecutive speeds. 

A surprising aspect of ITheorem 1.71 is that there is a positive probability (1/6) that Uq = Ui, even 
though each is uniform on [-1, 1]. Indeed, for any two particles there is a positive probability that their 
speeds are equal. This phenomenon can be thought of as a spontaneous formation of "convoys": sets of 
particle that have the same asymptotic speed, so their trajectories remain close. Our next result gives a 
full description of such a convoy. 

Theorem 1.8. Let the convoy ofO be Cq = {j : Uj = Uq}, i.e. the set of all j with the same speed as 0. 
Then Cq is p-a.s. infinite with density. Moreover, conditioned on Uq, Cq is a renewal process, and the 
non-negative elements o/Cq have the same law as the times of last increase of a random walk conditioned 
to remain positive, with step distribution 



The "times of last increase" of a walk Z are those indices n for which m > n implies Z,„ > Z„. In 
particular the convoys are infinite and they provide a translation invariant partition of the integers into 
infinitely many infinite sets. The convoys are essentially the process with density for second class 
particles, seen from a second class particle, as studied by Ferrari, Fontes and Kohayakawa in [5]. 



P(X = 1) = P(X = -1) = 




P(X = 0) = 




4 



2 



The TASEP speed process 



Amir, Angel & Valko 



6 



1.3 The ASEP 

As the name suggests, the totally asymmetric simple exclusion process is an extremal case of the asym- 
metric simple exclusion process: the ASEP. The ASEP is defined in terms of a parameter p e (1/2, 1], 
with p = I being the TASEP While most quantities involved depend on p, the dependence will usually 
be implicit. 

In the ASEP particles jump one site to the right at rate p e (I /2, 1] and to the left at rate p = I - p 
(we use the convention x = I - x). The generator of this Markov process is 

Lm = pificTnn) - m) + p{f«v) - m) o) 

n 

where cr„ and cr* sort the values in n, n + 1 in decreasing and increasing order respectively. 

While some of the questions above make sense also in this setting, there is a key difficulty in that the 
analogue of [Theorem 1.1 1 for the ASEP — conjectured below — is still unproved. Using the methods of 
Ferrari and Kipnis [7] it can be proved that Xo(t)/t converges in distribution to a random variable uniform 
in [-p,p], where hereafter we denote p = 2p - I. Note that the particles in the exclusion process try to 
perform a random walk with drift p (and they cannot go faster than that), that explains why the support of 
the limiting random variable is changed. In fact, in many ways the ASEP behaves similarly to the TASEP 
slowed down by a factor of p. 

Conjecture 1.9. In the ASEP, lim;^oo^o(0/^ exists a.s. (and the limit is uniform on [-p,p]). 

By the discussion preceding Corollary 1 1.2 1 this is equivalent to the following: 

Conjecture 1.10. The ASEP speed process measure p^^^'' is well defined and translation invariant with 
each Un uniform on [-p,p]. 

In order for statements about ihtASEP speed process to make sense we must assume this conjecture, 



and so some of our theorems are conditional on Conjecture 1.9 It should be noted that with minor 



modifications our results also hold assuming a weaker assumption, namely a joint limit in distribution of 
{X„(0/?}„ez- In that case, the speed process measure is still defined, even though the particles may not 
actually have an asymptotic speed. 

As noted there, [Theorem 1.41 holds also for the ASEP, with no additional condition. [Theorem 1.51 
becomes conditional: 



Theorem 1.11. Assume Conjecture 1.9 holds. Then p^ is a stationary measure for the ASEP: the 



unique ergodic stationary measure which has marginals uniform on [-p,p]. 

As in the case of the TASEP, this can be interpreted as follows: If an ASEP is started with initial 
configuration in [-p,p]^ with distribution p^^^^, then at any time t > the distribution of the process 
is also given by p^^^^ . Note that both the dynamics and the measure p'^^^'' depend implicitly on the 
asymmetry parameter p. 

A useful tool in studying the speed process is the understanding of the stationary measures of the 
A^^ + 1 type TASEP in terms of a multi-line process described below, developed by Angel [1] and Ferrari 
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and Martin [8]. There is no known analogue for these results that describes the stationary measure of 
the multi-type ASEP. Thus we need to use other (and weaker) techniques to extract information about 
the marginals of the ASEP speed process. This explains the contrast in the level of detail between the 
following results and the corresponding theorems above about the TASEP. 

Theorem 1.12. We have the following limit: 

\imV(Xo(t)<X,(t)) = ^^. 

Theorem 2.3 of [6] proves that the probability that particles and 1 interact at least once (i.e. one of 
them tries to jump onto the other) is In the next section we will show that this is equivalent to the 
just stated theorem. 

Our next theorem provides information about the joint distribution of {Uo, Ui}, assuming Conjec- 
ture [TH holds. 



Theorem 1.13. Assume Conjecture 1.9 holds. Let the measure i^'^^^ on [-p,p]^ be the marginal of {Uq, Ui} 
under fx^^^^. Denote by fi'-^^ the reflection off/^^ about the line x = y. Then on {{x,y) : —p<x<y<p] 
we have _ 

p.^i2,-p.^i2>^y^dxdy. 

We finish this section with a statement concerning the case Uq = Ui. Consider the total amount Jjj 

P l(\Xi(t) - Xj(t)\ = I) dt 

Theorem 1.14. In the TASEP, 7o,i = if and only ifUo = Ui. If Conjecture \1.9\ holds then the same 
holds for the ASEP. 

In the TASEP 7o,i = °° implies that there is at least one interaction between and 1 which means 
that they are a.s. swapped. (See the next section for a more detailed discussion.) Thus if Uq = Ui then 
eventually Xo{t) > Xi{t). In fact, this holds for any two particles in the same convoy: in ILemma 9.91 
we will prove that in the TASEP, particle will eventually overtake all the particles in its convoy with 
positive index. 



1.4 Overview of the paper 

The rest of the paper is organized as follows. [Section 21 provides some of the background: construc- 
tions of the processes and the multi-line description of the stationary measure for the multi-type TASEP. 
[Section 31 includes the proof of the symmetry property (ITheorem 1.41) and [Section 41 proves the station- 
arity of the speed process (Theorems 11.51 and 11.111) . Sections [6] and |7] include the results about various 
finite dimensional marginals of the TASEP speed process. [Section 81 deals with the proof of [Theorem 1.8[ 
Finally, in [Section~9| we prove our results on the ASEP speed process. 
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2 Preliminaries 



2.1 Construction of the process 

There are several formal constructions of the TASEP and ASEP. The one that best suits our needs seems 
to be Harris's approach [12]. We include the construction since there are several variations and the exact 
details are used in some of our proofs. The process is a function Y defined on Z x 1.+. Y^it) will denote 
the class of the particle at position k at time t. The configuration at time t is Y(t) = {Yk(t)]kez- The classes 
of particles will be real numbers, hence the configuration at any given time is in R^. Setting t = gives 
the initial configuration 7(0). 

We define the transposition operator t„, acting on by exchanging y„ and Yn+i, while keeping all 
other classes equal. Using this we can alternately describe the sorting operator cr„ by 



o-nY = 



1 Y otherwise. 



Thus cr„ has the effect of sorting Yn, Y„+i in decreasing order, keeping other classes the same. 

The TASEP is defined using the initial configuration and the location of "jump" points. The prob- 
ability space contains a standard Poisson process on Z x R"^, i.e. a collection of independent standard 
Poisson processes on R+, denoted r„. If {n,t) is a point of r„, then at time t the values of Y„{r) and 
y„+i(r) may be switched. In the TASEP they are sorted, i.e. Y(t) = Y{r) ■ cr„. This can be described as 
applying each of the operators cr„ at rate 1 independently. A simple percolation argument shows that this 
dynamics is a.s. well defined. (For any fixed t > there are a.s. infinitely many integers n so that there are 
no Poisson points on {n} x [0, t] which means that to define the process up to time t it suflaces to consider 
finite lattices.) 



The ASEP. Defining the partially asymmetric exclusion process requires additional randomness. Given 
the parameter p e (1/2, 1], we attach to each point (n, t) in the Poisson process an independent Bernoulli 
random variable , with P(X„ f = 1) = p. We can now define the probabilistic sorting operator p„ as 
follows: 

y^io-nY ifX„, = l 

\<y ifx„, = i. 

Thus with probability p the smaller classed particle is moved to the right position and with probability 
1 - p it is moved to the left position. When such an event happens we say that Yn{t) and F„+i(0 have an 
interaction (regardless of whether they were actually swapped). Note that if particles j interact in this 
way, then their order after the swap is independent of the order before the swap. The key observation is 
that after / < j interact in this way at least once, / has probability p of being to the right of j, and this is 

g l{\Xi{s)-Xj{s)\ = \)ds{ihe. 

total time /, j spend next to each other until time t) then 



P(XKO < Xj{t)\Ji^j{t)) = e'^^'i^^ + p{\ - e-^'-^^'^) = p + pe-'^'-J^'^ 



(4) 
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where the expression on the right is just the probability that there were no interaction between i and j 
until time t plus the probability that there was some interaction and at time t particle i is to the left of j. 
One of the consequences of @ is that 

lim F(Xi(t) < Xj(t)) =p+ p-Ee'^-'-' . (5) 

Thus [Theorem l.lll implies p + pEe~^ = ^ which gives 1 - Ee"-'" ' = But 1 - Ee"-'" ' is exactly the 
probability that there is at least one interaction between and 1 which shows why Theorem 2.3 of [6] 
and our [Theorem 1.1 21 are equivalent. 

In the TASEP case if there is an interaction between / < j then Xi{t) > Xj(t) after that. Thus in that 
case from ^ we get 

P(eventually > Xjit) | 7,- y = oo) = 1 

which explains the remark after [Theorem 1.141 

There is an alternate construction for the ASEP, which will be used in Section [31 Consider a Poisson 
process with lower intensity p on Z x R."^, but whenever it has a point (n, t) we apply at time t the operator 
nn rather then p„, where n„ is defined by 

TnY Yn > Yn+i with prob. ^ = (1 - p)lp 
Y Yn> Yn+i with prob. q = (2p - 1)1 p. 

Thus if the pair is in increasing order it is always swapped, while if it is in decreasing order it is swapped 
only with probability q. It is easy to see that every possible swap occurs at the same rate in the two 
constructions, hence the resulting processes have the same generator. 

2.2 Stationary measures for the multi-type TASEP 

The following theorem can be proved by standard coupling methods (see e.g. [11] where the same theo- 
rem is proved for the 2-type TASEP). 

Theorem 2.1. Fix every < Ai,...,An < 1 with J^Ai = 1. There is a unique ergodic stationary 
distribution vx for the N-type TASEP with 'P{Yq = k) = At- The measures Vx are the extremal stationary 
translation invariant measures. They are the only stationary translation invariant measures with the 
property that for each k, the distribution of {1[Y„ < k]}„^z is product Bernoulli measure with density 

2 j<k ^j- 

For the ordinary TASEP (with particles and holes) this stationary distribution is just the product 
Bernoulli with a fixed density. If we have an A'^ + 1-type TASEP then the structure of the stationary 
distribution is more complicated. The first description of vx for N = 2 was given by the matrix method [3]. 
[5] gave probabilistic interpretations and proofs of the measure and its properties. Recently combinatorial 
descriptions of vx have appeared as well. The (2 + l)-type TASEP was treated by Angel [1] (and see also 
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Duchi-SchaefFer [4]). These results were extended for all A'^ by Ferrari and Martin [8]. They give an 
elegant construction of v^i using systems of queues. 

We will now briefly describe the A^-line description of for the (A^ + l)-type TASEP. The 2-line 
case suffices for most of our results, with the exception of the results of Section U\ For a more detailed 
description and proofs see [8]. 

From here on we shall fix the parameters Ai, . . ., /Iw+i. Consider N independent Bernoulli processes 
on Z denoted 61,62, . . . , 6Ar where Bt has parameter Yji<k (these are the lines). From these lines we 
construct a system of A'^ - 1 coupled queues. The lines give the service time of the queues, and the 
departures from each queue are the arrivals to the next queue. 

It is important to observe that the time for the queues goes from right to left, i.e. 6,(n) is followed by 
6,(n - 1) and so on. The resulting system of queues is positively recurrent, so it can be defined starting 
at 00 and going over the lines towards -00. 

The z'th queue will consist of the particles that departed from the z'th line and are waiting for a service 
in 6,+i. This queue will consist of particles of classes {1, . . . , z}. When a service is available in 6,+i the 
lowest classed particle in the z'th queue is served and departs (to the next queue). If the queue is empty 
then a particle of class z + 1 is said to depart the queue. The departure process of each queue (i.e. the 
times and sequence of classes of departing customers) is the arrival process for the next queue. 

It is convenient to think of an additional queue with 61 as its service times. This queue has no arrivals 
(so it is always empty). The unused services introduce first class particles, which join the second queue 
whenever there is a service in 61 . These operations are evaluated for each n from line 1 to line in order. 
Let Qij(n) be the number of particles of type j in the z'th queue after column n of the multi-line process 
has been used. 

Note that each queue has a higher rate of service than of arrivals, so the queues sizes are tight, and 
the state with all queues empty is positively recurrent. In practice, the z'th queue has z types of particles in 
it, so the whole system of queues is described by ^Mzil non-negative integers. 

Theorem 2.2 (Ferrari-Martin), va is the distribution of the departure process o/Bn, with class N + \ (or 
empty sites) at those n when there is no service. 

As an example, and to clarify the graphic representation we use later, consider the following segment 
of a configuration of the 3-line process for n = {1,2,3,4}. Suppose both queues are empty at time 5. 
(This is denoted by the 0, exponent.) Here, O denotes a in the corresponding line, and • a 1 . Later, in 
cases where we do not care about a specific value we may use ® to denote that. 



At time 4, reading the rightmost column from top to bottom, there is no service in 61, so no first class 
particle joins the second queue, which therefore remains empty. There is a service in 62, and no particles 
in the first queue, so a second class particle joins the second queue. There is service in 63, so the second 
class particle departs immediately. Thus at time 4 the queue states are (0, 0). 

At time 3 a first class particle arrives to the first queue, and stays there since there is no service in 
the second queue. There is no further service in column 3, so the state at time 3 is ({1}, 0). There is no 
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departure, which is denoted by a 4 (or hole). At time 2 another first class particle arrives, and there is 
no particle in the second queue so the service in B3 gives rise to a third class particle departing. The 
states are now ({1,1}, 0). Finally, at time 1 a first class particle is served at both B2 and B3, departing and 
leaving queue states ({1}, 0). The resulting segment of is (1, 3, 4, 2). 



3 Symmetry 

Recall the operators ;r„ defined above. These act randomly on configurations, and the ASEP can be 
defined by applying each of the Markov operators n„ at rate p. 

Formally, nt is defined as acting on A1(5oo): probability measures on 5oo- Given a measure v on 5 00, 
we let TTnV be the distribution of n„ applied to a sample from v. Since r, and cr, also act naturally on the 
measures (in the same way), one finds the operator relation 

Tii = qTi + "qcTi. 

Note that p = I gives = so in that case tt, = cr,. In the case p = 1/2 we get ^ = 1 and ;r, = r,, so the 
process reduces to a symmetric random walk on 5oo- 

The crucial observation leading to [Theorem 1.41 is the following lemma. 

Lemma 3.1. Fix any p > \/2, and sequence i\, . . in- Then 

ni„---7Ti^ - id = (;r,-, •••;r;„ •id)"\ (6) 

That is, applying a sequence of ;7r, 's in the reverse order to the identity leads to the inverse permutation. 
This is trivially true when p = \/2 and n = t, but requires proof for other p. When p e {1/2, 1} the 
operator is deterministic and this distributional identity is an equality of permutations. 

Proof of lTheorem T~4\ The theorem follows from ILemma 3TT] since at any finite time at each i there is 
positive probability (e~') that no swap has occurred. Each such i separates Z into two parts with inde- 
pendent behaviour, so the state of the process is a product of finite, mutually commuting permutations. 
The distribution of the sequence of applied operators between such inactive locations is symmetric in 
time. □ 

We now prove [Lemma 3 . 1 [ In the case of the TASEP. ILemma 3. H and [Theorem 1.4l were first proved 
in [2]. To prove the lemma in the general case, we start with the following facts about the transposition 
operators. The identities are readily verified, and the last claim is known as Matsumoto's Lemma (see 
e.g. [?, Theorem 3.3.1]). 

Fact 3.2. The operators ti satisfy the relations 

= / (7) 

'^i'^j = '^Fi for \i - j1 > 1 (8) 

TiT^Ti = T;+lT,T;+i (9) 
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where I denotes the identity operator. With these relations the operators {r,} generate the symmetric 
group. Furthermore, it is possible to pass between any two minimal words of the same permutation using 
only relations dS]),®. 

The ;r's satisfy similar relations: 

Lemma 3.3. The operators {ni] satisfy the relations 

= ql + qm (10) 
niTTj = njnj for\i- j\> I (11) 

1 = i^i^i+i- (1 2) 

Note that only the first relation differs from the corresponding relation for t. When p = 1/2 these 
reduce to the relations for r. In the case p = I the first relation becomes cr^ = cr,. In that case, the only 
non-trivial relation is (fT2l) which is true since both sides have the effect of sorting the three terms involved 
in decreasing order. 

Proof. (flOl) is easy to check, and (fTTI) is trivial. For (fT2l) . using n = qr + qcr and expanding, we need to 
show that 

<?^(t,T,+iT;) + q^q(TiTi+iO-i + TiO-i+iTt + Cr,T,.+ iT,) + q'qi^TiO-MO-i + Cr;T,+iCr; + Cr,Cr;+iT,) + 'q{(Ti(TM(T^ 

is unchanged by exchanging / and / + 1 . It is easy to verify that 

TiT,+iT, = T,+iT,-T,-+i r,-T,-+icr = (T ,+iT,T,-+i r,-cr ;+ir,- = T,-+icr ;r,-+i 

o-iO-i+ 1 o-i = ai+ 1 cr,cr,+ 1 cr;T,+ 1 r,- = r,+ 1 T,cr,+ 1 

so it remains to show 

T,cr,+i CT; + cr,r,+icr,- + cr,cr,+iT,- = r,+iCr,cr,.+i + cr,+iT,a-;+i + cr,+iCr,T,+i 

We may assume / = 0. Since only the relative order of ?7o, ?7i, '72 matters, we may assume these are {0, 1, 2} 
in some order. Apply these operators to the 6 orders to get the following table. In each column, the entries 
in the top half are a permutation of the entries in the bottom half, so adding the first three operators is the 
same as adding the last three. 





012 021 102 120 201 210 


TqO-iCTq ■ 7] 
O-qO-iTq ■ 7] 
O-qTiO-q ■ 7] 


210 120 210 120 120 120 
210 210 201 210 201 210 
210 210 210 201 210 201 


TiO-qO-i ■ 7] 
O'iO-qTi • 7] 
O-iToCTi • 7] 


210 210 201 201 201 201 
210 120 210 120 210 210 
210 210 210 210 120 120 



□ 
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Proof of VLemma 3.1\ Let X = t,, • • -r,,,. If this is a minimal word for X in ^oo, then applying the nis 
results in X with probability 1. Since the reverse word is a minimal representation for the lemma 
holds in this case. 

Consider the set of words that satisfy the lemma. It contains all minimal words. The proof now 
proceeds by induction. Take some non-minimal word, and let k be maximal such that X = t,j • • • t,-;, 
is minimal word, let X = tj^ ■ ■ -Tji^ be another representation of X with jk = ik+i- By the induction 
hypothesis, the claim holds for (ji, . . . , jk-iJk+i, ■ ■ ■, h) and (ji, . . . , jk-iJk+i, ■ ■ - , in)- By relation 1 we 
can replace n^^ by a linear combination of ;7r,j and /, so the claim also holds for (ji, . . . , jk, ik+i, ■ ■ ■ , in)- By 
applying relations (fTTI).(fT2l) to the first k terms we find that our original word also satisfies the claim. □ 

Note: The proof actually shows that any word in the ;r's can be reduced (as an operator) to some 
convex combination of words corresponding to minimal words. 

Corollary 3.4. Consider the infinite type TASEP with initial condition y„(0) = n. Then {-j^}„gz converges 
weakly to n as t ^ oo. 

Proof. For any t this process has the same law as {-j^}„6z, which converges a.s. to a process with law 
fi. □ 



4 Stationarity 

We will give two different proofs of the stationarity of the distribution of the speed process. The first 
is specific to the TASEP, and is reminiscent of coupling from the past. It uses the Harris construction 
directly. The second proof is based on the symmetry between {X„{t)} and {Yn(t)} (or more specifically 
Corollary 13.41) . The second proof holds also for the ASEP, word by word, under the assumption that 
Corollary 13. 41 is true for the ASEP (which is weaker then Conjecture |1.9l ). 

4.1 Coupling proof 

Lemma 4.1. Consider two TASEPs Y, Y' defined via the Harris construction as the function of the same 
Poisson process onZx M^. We set the initial conditions as Y^iO) = n and Y'{0) = (TqY{0) (i.e. particles 
and 1 are switched initially in Y'). Let {Un} = {\imt^ooX„(t)/t} denote the speed process of Y, and 
{U'J = {limt^oo X'^(t) / 1} denote the speed process ofY'. Then U' = ctqU. 

Proof. All particles other than {0, 1 } are either larger or smaller than both and 1 , so any swaps involving 
a particle other than {0, 1} will occur or not occur equally in Y and Y' . It follows that for any i i {0, 1} 
we have Xj(f) = X.(t)and hence Ui = U-. Similarly, since and 1 must fill the only vacant trajectories, 
{Uo, Ui) = {U'q, U[} as an unordered pair. 

In Y' particle is always to the right of particle 1, so U'^ = max{Uo, Ui} and U[ = mm{Uo, Ui}, 
completing the proof. □ 



The TASEP speed process 



Amir, Angel & Valko 



14 



Proof of fTheorem L5\ using coupling. Consider a Poisson process on Z x R. Half of the process, namely 
the restriction to Z x is used in the Harris construction of the TASEP. Similarly, for any s eRwe can 
translate the Poisson process by s (i.e. take all points of the form (n, t + s) where (n, t) is in the original 
process), and take the restriction to Z x R"^, which can be used in the Harris construction to get a different 
(though highly dependent) instance of the TASEP. 

Let Un{s) be the speed process resulting from the Harris construction using the translated Poisson 
process. Clearly for every s, U(s) has the same law /i, so we are done if we show that Un(s) evolves as 
a TASEP (with time parameter s). Consider the effect of an infinitesimal positive shift s. The shift adds 
new cr operations, to be applied before the original sequence of operations. These are added at rate 1 at 
each location. By the previous lemma, the effect on the resulting speeds of applying cr„ before using the 
same Poisson process is to apply cr„ to the speeds, which is exactly what we need. □ 

It is interesting to note that in the Poisson process Z x R, the part on Z x R^ is used to determine the 
"initial" speed process U(0), and the restriction to Z x R" is used exactly as in the Harris construction to 
generate the TASEP dynamics of U{s). 



4.2 Symmetry based proof 

Proof of Theorems \1.5\ and \l.ll\ using symmetry. We write the proof for jjL, but it holds verbatim for /u^^^^ 



under Conjecture 1 .9 



Let Ps be the evolution operator for the Markov process corresponding to the generator L on R (see 
([U)). To prove stationarity it is enough to show that for every < s and every bounded continuous local 
function / : R^ ^ R we have 

J PJdfi = J fdiu. (13) 

Consider the process {y„(0}«ez started from Y„(0) = n and denote the distribution of \ ^^\ by v,. By 
Corollary 13.41 the weak limit of v, is fj. which means that for every local bounded continuous function 
/ : R^ ^ R we have 

By the continuity of P^ we get that for any fixed s 



f PJdv, > f PJdfX. 

J '^'^ J 



But J Psfdvt = '-Y j fdvns which (for a fixed s) converges to J fdji. Thus (fT3] ) and the theorem follow. 

□ 



5 Basic properties of stationary distributions 

In this section we present a medley of simple results concerning the (T)ASEP and its stationary distribu- 
tions. These are only weakly related to each other, and are collected here for convenience. 
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Proposition 5.1. fi is ergodicfor the shift. Under\Conjecture 1.9\so is fi 



,ASEP 



Proof. Consider the setup of Corollary 11.21 and use the Harris construction with independent standard 
Poisson processes T„ on the interval [0, oo) to define and the variables Xn{t). Then the limit process 
{Un}nez is measurable with respect to the cr-algebra T generated by the i.i.d. processes T„ (n e Z). Since 
'F is generated by i.i.d. processes any translation invariant event in has to be trivial. But then the same 
thing must be true for any translation invariant event in the cr-algebra generated by {U„}nez as this is a 
sub-cr-algebra of !F. □ 

There are three possible "reflections" for the ASEP. One may reverse the direction of space, so that 
(low classed) particles flow to the left and not right; one can consider the time reversal of the dynamics, 
and one can reverse the order of classes, (or keep the same generator but replace class k with -k, or 
N+l-k, etc.) It is easy to see that reversal of both space and class order preserves the original dynamics. 
This is called the space-class symmetry of the TASEP/ASEP. 

The following proposition is the space-class symmetry of the speed process, and follows directly 
from the corresponding symmetry of the ASEP process. 

d 



Proposition 5.2. For the TASEP {U„}„^z = {-U-„}nez- This also holds for the ASEP, assuming \Conjecture 1.9 
holds. 

The following observation and its corollary provide an important connection between the distribution 
of the speed process and the stationary measures of multi-type ASEP. These connections will be used to 
extract information on the joint distribution of the speeds of several particles in Sections [6] and |71 

Lemma 5.3. Let {rin{t))nj be an ASEP, and let F : R ^ R be a non-decreasing function. Then {F{T]„(t))]„^, 
is also an ASEP (with the same asymmetry parameter). 

Proof. The ASEP is defined as applying to ri(t) each of the operators independently at rate 1 . Applying 
a non-decreasing function to each coordinate commutes with every ni, hence {F(T]„(t))}n t is just the ASEP 
with initial configuration {F(?7„(0)}„. □ 

Corollary 5.4. If F : [-1,1] {l,...,N} is non-decreasing then for the TASEP the distribution of 
{F(Un)} is the unique ergodic stationary measure of the multi-type TASEP with types {1,...,A^} and 
densities Ai = j fieb(F"^(/)). 



This also holds for the ASEP (and its corresponding multi-type stationary measure) under Conjecture 1.9 



Proof Let pf denote the distribution of {F{Un)}. Since p is ergodic, so is pf- The marginals are as 
claimed since each t/„ is uniform on [-1, 1]. 

To prove that pp is stationary, start a TASEP Y„{t) with initial configuration y„(0) = U„. By lLemma 5.31 

{F{Y„{t))}„^t is a N-type TASEP. Since p is stationary, Y(t) also has law p, and so {F(Yn(t))}nj = {F(Y„{0))}„ t, 
hence pf is also stationary. 

The result for the ASEP follows the same way. □ 
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The next proposition shows that a TASEP started with uniform i.i.d. classes must converge to the 
speed process. In particular, even though classes in the i.i.d. initial distribution are a.s. all different, the 
process converges to the speed process which has infinite convoys of particles with the same class (see 
Section [8]). Thus the TASEP dynamics has the effect of aggregating particles with increasingly closer 
speeds next to each other. 

Proposition 5.5. Consider a TASEP where Y„{Qi) are i.i.d. uniform on [-1, 1]. Then {F„(0},iez converges 



weakly to fx. The same holds for the ASEP under Conjecture 1.9 



Proof. Let Vt be the distribution of Y{t) for the process Y of the lemma. We need to show that f gdvt > 

f— »oo 

J gd^ for any fixed bounded and continuous function ^ : [-1, 1]^ ^ R.. 

If we start the A^-type TASEP with an i.i.d. product measure initial distribution then its distribution 
converges to an ergodic stationary measure with the same 1 -dimensional marginal. (This can be shown 
by standard coupling arguments introduced by Liggett, see e.g. [11] or [12, Chapter 8].) 

Using Lemma [53] and Corollary 1 1.61 it follows that for any non-decreasing step function F on [-1, 1] 
the process {F(y„(0)}« converges in distribution to {F(i7„)}„. 

For an integer M let Fm{x) = which maps [-1, 1] to {i/M, i e [-M, M - 1]}. Define the operator 
F^ on configurations, as the operator that applies Fm to each coordinate: F^(r])n = FMiUn)- Since g is 
continuous we can select M such that ||^ - ^ o < s. By the triangle inequality we have 

J gdv,- J gdiu\<2s+\J goFldv,- J goF^dfi 

g o F^ is g applied to a TASEP with finitely many types, so it can be made smaller than s by taking t 
large enough. □ 



6 Two dimensional marginals of the TASEP speed process 



The key tool for analyzing the joint densities of the speed process is Corollary 5.4 This states that if the 
speed process is monotonously projected into { 1 , . . . , ^, ^ -i- 1 }, then the result is the stationary measure of 
the multi-type TASEP with suitable densities. In the TASEP, the latter is given in terms of the multi-line 
process (see Subsection 12.21) . More explicitly, we will use the following projections, to which we refer 
as canonical projections. Let x = (xi, . . . ,Xk)he an increasing sequence taking values in [0, 1], with the 
conventions that Xq = and Xk+i = 1. Define F : [-1, 1] {I, . . . ,k,k + 1} by 

\ + u 



F{u) = FJja) = min{z : u < Xi], where u 

L 

Note that if u is uniform on [-1, 1] then F(u) = i with probability Xi - x\-\. Let V, = F{Ui), so each V,- 
has distribution controlled by the x's. It is not hard to see that the cr- field generated hy V\, . . . ,Vk (or any 
k fixed indices) for all possible x's is the same as the cr-field oi Ui, . . . ,Uk- 

The scheme of our argument should now be clear. The distribution of V is given by a multi-line 
process, and can be computed explicitly. Considering the resulting probabilities as functions of x allows 
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us to recover the joint density of the corresponding speeds. This last step is done by taking suitable 
derivatives w.r.t. jc,'s to get the density. In order to find the joint density of k particles we work with the 
k-line process. In this section we use this approach to prove results about two dimensional marginals 
of fi. We prove [Theorem 1.7] which gives the joint distribution of {Uq, U\) and generalize this result for 
the joint distribution of any two speeds. In the next section we give some results for higher dimensional 
marginals. 



6.1 Two consecutive speeds: Uq,U\ 

Proof of lTheorem 7771 We compute the probability that Vi = 2 and Vq is each of 1, 2, 3 (recall that as the 
highest class particles, 3's are equivalent to holes). The queue of the two line process is a single, simple 
queue, so indices are not needed. In order to have a second class particle at position 1 we need an unused 
service. This means the queue must be empty: Q{2) = 0, and there must be a particle at the bottom line 
but not at the top line in position 1 . The intersection of these events has probability X2 - xi (as this is 
the density of second class particles). More importantly, they depend only on the 2-line configuration in 
positions {1, . . . , oo}. Since on this event the queue is also empty at position 1, the class Vq depends only 
on the 2-line configuration at position 0. 

In particular, to get a first class particle, Vq = I, the only possibility is to also have particles in both 
lines in position 0. This leads to 



^(Vo = 1, Vi = 2) = P(||f ) = X,X2(X2 - 



Xi). 



We shall also denote this probability by yUx(1.2) for compactness, as this is the probability of seeing 
consecutive particles of classes 1, 2 in the stationary measure //x- Similarly we have 

fj.^{2, 2) = P (SS^) = ^1 ^(-^2 - -^i), 
fi.(3,2) = ¥(^^=Mx2-x,). 

Here, ® indicates no restriction on the top line in that position and y = I -y. 
To calculate the densities of the two speeds we find for example 

P(t/0 < 2jCi - 1 < t/i < 2X2 - 1) = A'x(l, 2) = XiX2(X2 - Xi). 

Thus to find the density at {uq,U\) for uq < Ui we need to take derivatives w.r.t. X2 and xi, and set 
^2 = (1 -I- Mi)/2, xi = (I + Mo)/2. Remembering the Jacobians (1/2) we find 

(1 \ / 1 \ ^1 — ^0 

jdxi 1^2^x2) l^x(i, 2) = — - — duodui for Mo < mi. 

Similarly, to find the density at (mq, mi) for mq > mi noting that the Jacobians now have reversed signs 
we find ^ 

F(Uo 6 duo, Ui 6 dui) = {-jdxi)i^-jdx2)l^x(^,2) = -duodui for mq > mi. 
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Finally, to find the (singular) density along the diagonal, consider /ix(2, 2) and let X2, xi ^ We 
have 

P([/o, C/i edu) = \ lim ^""^ ' = ——du. □ 



6.2 Two distant speeds: Uq, Uk 

The two line process also yields formulae for the joint density of two distant particles. However, the 
result is not as compact as for the case of two consecutive particles. 

Theorem 6.1. For any k > Owe have: 

• The joint density ofUo, Uk on {Uq > Uk) is 1/4 (so ¥{Uq > Uu) = 1/2). 

• On {Uq < Uk} the density is a polynomial of degree 2k - 1. 

• On the diagonal {Uq = Uk} the density is a polynomial of degree 2k. As k ^ oo, the density on the 
diagonal {(w, m) : |m| < 1} is asymptotically yjj^- 

It is possible to prove exponential convergence of the density on {Uq < Uk} to 1/4, though we do not 
pursue that direction here. The fact that as ^ ^ oo the distributions of Uq and U^ become independent 
follows from ergodicity, or can be read from (fT6l) below. 

The theorem follows easily from the next two lemmas. Let {S,,} be a random walk with sites in {±1,0} 
with probabilities {p+, pq}, and denote M„ — max,<^ lb 

Lemma 6.2. Fix < x < y < 1, and let S,,, M„ be as above with 

p+ = xy p- = xy pq = xy+ xy. 

Then we have the following: 

¥(x<Uu<y<Uo) = {y-x)y, (14) 
P(i7o, Uu 6 {x,y]) = {y- x)xyP(M,.i = 0) (15) 
¥(U(i<x<Uk<y) = {y- x)xy + {y-x) xy¥{Mk-\ > 0) (16) 

Note that the steps of S are the difference of two Bernoulli random variables, and therefore S j = 

prob. 

Bin(7, .x:) - Bm( j,y). In particular, for any fixed x < y v/e have S j > -oo, and asymptotically the 

speeds are independent. 



Proof. By Corollary 5.4[ ¥(x < Uk < y < Uq) = /i.v,v(3,2) (where p^-^y the extremal stationary 3 type 



TASEP with densities x,y - x,\ - y). Using the 2-line description of fXx,y we have Vq = 1, V^t = 2 if and 
only if we see the two line configuration 
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Having the hole in the bottom line at position has probability y and this is independent of having a 
second class particle at position k. 

Similarly, to have Uq, 0^ & [x,y] we need the configuration 



o- • 


■o 


••■ 


■• 



with intermediate configuration leaving the queue empty at position 1. Let S j be the number of particles 
in the top line in positions {!,...,]} minus the number of particles in the bottom line in those positions. 
The condition that the queue ends up non empty is equivalent to {maxo<y</t-i S j > 1}. The claim follows. 

Finally, the third case follows from the first two since the three probabilities must add up to P(f/i € 
[x,y])=y-x. □ 

Lemma 6.3. LetS„, M„ be as above with p+ = p.. Then P (M„ = 0) = P(5„ 6 {0,-1}) 

Proof. Reflection at the hitting time of 1 shows that 

P(M„ > 0,5„ < 0) = P(M„ > 0,S„ > 2) = P(5„ > 2) = P(5„ < -2). 

It follows that 

P(M„ > 0) = P(5„ > 0) + P(M„ > 0, S„ < 0) = 1 - P(5„ 6 {0, -1}). □ 

Proof of lTheorem 6J\ The case Uq> U\ is just the double derivative of (fT4l) . 
For the case Uq = U\, note from (fT5l) that the density along the diagonal is 

x,y^u 2(y - x) 

where M^^^i is the maximum of a symmetric random walk with p+ = = xx. Using the prior lemma, 
since p+ = p- we get 

P(M,_i =0) = P(5,^i 6{0,-l}). 

This is clearly polynomial. Using the local central limit theorem, P(5yt-i = a) ~ ^^i-^ for any a G 
{0,-1}, and our claims follow. 

For the case Uq < Ui, taking derivatives of (fT6l) shows that the density is polynomial as claimed. □ 



7 Multiple speeds 

In this section we will prove some results about the joint distribution of more than two speeds. In princi- 
ple, any finite dimensional marginal of the distribution can be derived from lTheorem l.Sl along the same 
lines as used above for the joint distribution of Uq, Ui. This gives the joint distribution in terms of the 
stationary measure of the multiple queue system. Some aspects of the joint distribution have particularly 
nice formulae and we proceed to present some of these. 
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1 . The next subsection determines the probability that out of the first n particle a given one is the 
fastest. 

2. The following result shows that the speed of a fast particle is independent from those of adjacent 
particles it overtakes. More precisely: if c 6 [-1, 1] then conditioned on the event that Uo > c and 
Ui, . . . ,U„ < c, the random vector (t/i, . . . , U,,) and Uq are independent. 

3. Next, we show that on the event {Uq < U\ < ■•■ < there is a pairwise repulsion between the 
particles: the density function is given by n\ times a Vandermonde determinant. 

4. Finally, we give the full description of the joint distribution of {Uq, Ui, U2). Their distribution is 
absolutely continuous with respect to the Lebesgue measure on each of the 13 subsets of [-1, 1]^ 
corresponding to a given order of these speeds (these include the cases where two or three speeds 
might be equal). In [Theorem 7.71 we determine the densities on all of these subsets. 



7.1 The fastest particle 

As a first example, we compute the probability that particle i will be the rightmost of {1, . . . , n} for all 
t > to. This proves and generalizes a conjecture of Ferrari, Goncalves and Martin [6] that the probability 
of particle overtaking particles 1 through n is Note that this is not quite the same as saying that Ui 
is the maximal of {f/i , . . . , [/„}. Due to lLemma 9.9[ this event allows f/, = Uj for j > i but not for j < i. 

Theorem 7.1. For any n and any k e [l,n] 

In 



lim P {Xuit) = max{Xi(0, • • • , X„m 



(n + k- \){n + k) 



Lemma 7.2. Let X = Bin(m, p) and Y = Geom(^) be independent binomial and geometric random 
variables. Then 

P(y <X) = \-q{p + pqT. 

Proof. We have ¥{Y > X) = Zi ("l)p' p'"''q'^^ =q(p + pqT ■ □ 

Proof of lTheorem 7J\ Note that since the index of the rightmost particle (of {I, ... ,n}) is non-increasing 
in time, the event in the statement is equivalent to particle k being the rightmost for all t > to for some to. 
Bv lLemma 9.9[ which we prove in Section|9l particle i eventually passes particle j for i < j if and only if 
Ui > Uj. Thus k will eventually be the rightmost particle of particles {1, ...,«} if and only if Uk > Ui for 
1 < / < k and Uic > Ui for k < i < n. Call this event Ej,. 

As an intermediate step we will compute the probability that this happens and U/c e du for some 
u e [-1,1]. Integrating over u will give the theorem. Fix x = (xi,X2), < a:i < .^2 < 1 and consider the 
event Ej^^g, that for all z e [1, n] we have that 



Uie 



[0, Xi] i < k, 
[xuX2] i = k, 
[0, X2] i > k. 
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Thus Ek,^ says that up to the partition resulting from the vector x, the event Ej, holds. 

Projecting into the 2 + 1 type TASEP using £^/t,x is mapped to the of event of having fc- 1 first class 
particles followed by a second class particle, followed by n - k particles of either class (but no holes). 
This requires in positions l-n a configuration of the following form: 

• • •••• • • •! ' 

where the first hole in the top line is in position k, and the size / of the queue can be no greater than the 
number of holes in the top line in positions {k+ I, . . . ,n}. Since the number of holes in the rest of the top 
line has the binomial distribution Bin(n - k, Jj") and the queue state is an independent Geom(=^j, we 
find after simplifying that 

F(Ek,^) = ir^^p|Geom|^^j < Bin(n - /t, = 4'^ ~\A ' A^i^ ~2 

(noting that q + pq of the previous lemma simplifies to X1/X2). 
Taking a limit as .)C2, xi — > a we find 

¥(Ek, Uk 6 dy) = lim ^^^^ = y^+i'-^ (^(^^ + \ - k) - (n - k)y) dy. 

Finally, integrating overy e [0, 1] gives 

¥(E,)= f f^'^-\{n+\-k)-{n-k)y)dy = - f , ■ □ 

Jo {n + k-\){n + k) 

7.2 Independence when swapped 

The following result shows that the speed of a fast particle is independent of speeds of adjacent particles 
that it overtakes. 

Lemma 7.3. Fix c e [-1, 1] and a measurable set A c [-1, c]". Then we have 

fu(Uo>c\(Uu...,U„)eA)=iuiUo>c). 

Furthermore, conditioned on Uq > c and (Ui, . . . ,U„) e Awe have that Uq is uniform on [c, 1]. 

Proof. Since products of intervals span the cr-field, it suffices to prove the analogous statement for the 
M-type TASEP (in fact M = n + 1 is enough). Consider a TASEP measure /ix where holes have density 
1 -c, so that speeds greater than c correspond to holes. We need to show that for any classes /i , . . . , z„ < M 

MVo = M\V,=iu...,Vn = in) = fi.iVo = M). (17) 

To show this we consider the multi-line process. There the classes of Vi, . . . , V,, are determined by 
the lines in positions [1,00). On the other hand, Vq = M requires only that Bm{0) = 0, hence the 
independence. 

To get the second claim, note that /i(t/o > = ^ th^t (flTI) also applies (with the same set A) for 
any c' > c. □ 
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Corollary 7.4. The (Ui,..., Un)-marginal offi has a constant density function 2 " on the set {Ui > • • • > 
U„}. 

Proof. The events that the speeds are in small intervals around the m;'s are independent. □ 
7.3 Repulsion when unswapped 

Here we derive the density function of the n + 1 -dimensional marginal of on the event {Uq < ••• < ?/„}. 
The result is given in terms of a Vandermonde determinant defined by 

a<i< j<b 

We start with a simple lemma about these determinants. 
Lemma 7.5. Let xq < • • • < x„. Then 



J J J Xi-i<yi<Xi 



Proof. We use the fact that A(y) is the determinant of the Vandermonde matrix: l^i^Jj) = det^jj^ _ . 
Since the determinant is is linear in the rows and each appears in a single row, we can integrate row by 
row to find 



rrr ^imwu = rrr n-=i = doJ^"'-' 

J J J Xi-i<yi<Xi J J J Xi-i<yi<Xi 

where M = (x'. , - x^T . Extend M to an (n + 1) x (n + 1) matrix M' by 

\ • i '/i,j=l 



= — detM, 
n\ 



1 


n ^ 

Xo • • • Xq 









M 








M' = 



Clearly detM = detM'. However, by sequentially adding each row to the one below it we find detM' 



x, = Ao„(x), completing the proof. 

' / i,j=0 ' 



Lemma 7.6. Let = xq < X\ < • ■ • < Xn < x„+i = 1, and be the corresponding n + I type TASEP 
stationary measure. Let Qn be the probability that all queues are empty at any specific location of the n 
line process. We have the following 



L p.^(2, ...,n)= fiQ/i G [2, n], Ui G x,]) = Ai,„(x), 
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2. yUx(l, ...,«)= nO/i e [1, «], Ui e x,]) = Ao,„(x), 

3. The density of U\, . . . ,Un on the event U\ < ■ ■ ■ < Unis n!Ai„(M). 



4. Qn = TTr. 



Ai.„(x) 



Proof. The proof is by induction on n. For n = I, claims 1 and 4 are trivially true, and 2,3 hold since the 
speeds are uniformly distributed. 

The key observation is that the only n-line configuration giving particles of classes 1, . . . , n is 



•00-- 


•oo 


••0-- 


•oo 




•oo 


•••• • 
•••• ■ 


••6 
••• 



(with all queues empty). Since the queue state is independent of the configuration in these positions, we 
find 



i=l 



This implies equivalence of claims 2 and 4. 

Similarly, the only configuration giving particles of types 2, . . . , n is 



00-- 


•oo 


•0-- 


•oo 




•oo 


••• • 
••• • 


••6 
••• 



This implies equivalence of claims 1 and 4 (since Ao,„(x) = Ai „(x) Y\ ^d- 

Next, we argue that claims 2 and 3 are equivalent. Claim 2 follows from 3 by ILemma 731 Claim 2 
also implies claim 3, since the density is the multiple derivative H/Li ^ of the probability of claim 2. 

Thus for any given n, the four claims are all equivalent. To complete the proof (by induction) we note 
that claim 3 for a given n implies claim 1 for n + I. This also follows from ILemma 7 . 5 1 in the same way 
as claim 2. □ 



7.4 Joint densities for 3 consecutive particles 

This section contains the complete description of the joint distribution of (Uq, Ui, Uj). The distribution 
is absolutely continuous with respect to the Lebesgue measure on each of the 13 subsets of [-1, 1]^ 
corresponding to a given order of these speeds (these include the cases where two or all three speeds 
might be equal). In [Theorem 7.71 we determine the densities on all of these subsets. 
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Theorem 7.7. The joint distribution of Uo,Ui,U2 is given by the following table, broken according to 
their order 



order 


density 


order 


density 


Mo < Ml < M2 


^(M2 - Mi)(Mi - Mo)(M2 - Mo) 


Mo = Ml < M2 


^(M2 - Mi)(l - Mj)(2 + 3m2 - Ml) 


Mo < M2 < Ml 


3^(M2 - Mo)(2 + 4mi - 3m2 - 3mo) 


Mo < Ml = M2 


^(Mi - Mo)(l - u])(2 - 3Mo + Ml) 


Ml < Mo < M2 


■^(M2 - Mo)(2 + 3m2 + 3mo - 4mi) 


Ml < Mo = M2 


^(M2 - Mi)(l - ul) 


Ml < M2 < Mo 


- Ml) 


Mo = M2 < Ml 


^(Mi - Mo)(l - M^)) 


M2 < Mo < Ml 


- "0) 


Ml = M2 < Mo 




M2 < Ml < Mo 




M2 < Mo = Ml 


1^(1 -"^) 


Mo = Ml = M2 









Proof. Fix < x\ < X2 < xt, < \. Define F = as above, and V, = F{Ui). To calculate the densities 
of the various simplices and facets, we calculate partly the distribution of V, and take suitable derivatives 
and limits. It is interesting to note that there are several possible class configurations for each case. For 
example the case {Uq < U\ < U2} can be deduced from each of /ix(l,2, 3), /Zx(l,2, 4), //x(l,3,4), and 
//x(2, 3, 4). Careful choice of the cases to consider can simplify the computations significantly. 

Not all cases need to be worked out. Space-class symmetry reduces many cases to others. [Theorem 1.7[ 
Lemmas 17.31 and l7!6l and Corollary 17 .41 imply several cases. Thus even though all 13 cases can be com- 
puted using this method, only 4 are essentially new. 

The following table summarizes the proofs for the 13 weak orders of Uo,Ui,U2. Here, A(x) = Ai 3(x). 



order 


V 


yUx(V) 


remarks 


Uo<Ui< U2 


1,2,3 


Xi.x;2.x:3A(x) 


Lemma 7.6 


Uo<U2< Ul 


2,4,3 


(Ji + ^)^A(x) 


new 


Ui<Uo< U2 


2,1,3 


Xi{X2 + Xi,)A(x) 


space-class symmetry 


Ui<U2< Uo 


4,2,3 


x^A(x) 


Theorem 1.7llLemma 7.3 


U2<Uo< Ul 


2,3,1 


XiA{x) 


space-class symmetry 


U2<Uy< Uo 


3,2,1 


Xl(X2 — Xi){X3 — X2) 


Corollary 7.4 


Uo = Ui< U2 


2,2,3 


'XiX2Xt,A(x) 


new 


Uo<Ui = U2 


2,3,3 


'XiX2X3A(x) 


space-class symmetry 


U:<Uo = U2 


3,2,3 


x^x^Aix) 


new 


Uo = U2< Ul 


2,3,2 


Y[x2A(x) 


space-class symmetry 


Ui = U2< Uo 


4,2,2 


X1X2 Xi,(X2 — Xi) 


Theorem 1.7llLemma 7.3 


U2<Uo = Ul 


3,3,1 


Xi X2Xt,{X3 — X2) 


space-class symmetry 


Uo = Ui = U2 


2,2,2 


Ji'xliX2-Xi) 


new:ITheorem 1.81 



The case {Uq < Ui < U2} is a special case of lLemma 7.61 while the case {U2 < Ui < Uq) is a special 
case of Corollary 7.4 The cases {Ui < U2 < Uq} and {Ui = U2 < Uq} follow from joint distribution 
of t/i, i72 [L7] together with lemma 1731 Each of the five cases {U2 < Uq < Ui}, {Ui < Uq < U2], 
{Uq < Ul = U2},{Uo = U2 < Ui}a.nd{U2 < Uq = Ui} follows by space-class symmetry ([Proposition 5.2[) 
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from the cases {Ui < U2 < Uo], {Uo < U2 < Ui}, {Uq = Ui < U2], {Ui < Uq = U2] and {f/i = U2 < Uq} 
respectively. 

It therefore remains to prove just 4 new cases: {Uq < U2 < Ui), {Uq = Ui < U2}, {Ui < Uq = U2} 
and{t/o = f/i = U2}. 

For the case {Uq < U2 < Ui},we compute //x(2, 4, 3). The only 3 line configurations that give these 
types are 



omo 



and 



000 

•®o 
mom 



f3 



Therefore 



/Zx(2, 4, 3) = Xi X2 X2x\ XT,{Xi+ X2) yUx(empty queues) 

_ — 2 —2 — .— ^—. ^13(X) 
— X\ X2 X2X^ Xt,\ Xi -r X2) 



Xi X2 X2X^ 



= X3 ( Xi + X2) Ai3(x) 
Taking derivatives we find the density of (/q, C/i, in the domain {Uq < U2 < f/i}is 

/Zx(2, 4, 3) = (mi - M2)(2 + 4mi - Suq - 3u2). 



-d 


-d 


-d 




8x3 


X3=ui 9X2 


JC2=«2 ^-^y 





A linear change of variables gives the formula in terms of u\,U2,uj,. 

For the case [Uq = U\ < U2}, we consider /ix(2, 2, 3). The only 3-line configuration giving this result 

is 



poo 

■ 'O 



10,0 



Thus 

yUx(2, 2, 3) = jCjjCj^xTVxCempty queues) = ir-'C2-^3Ai,3(x). 
Taking a derivative w.r.t. X3 and letting X2 X\ gives the density of the t/;'s to be 



km — 

x\,X2-^ua X2 — X\ OXj 



;Ux(2, 2, 3) = Uq Uq(U2 - Uq){3u2 - Uq). 



As above, a change of variables gives the claim. 

For the case {Ui < Uq = U2} we consider yUx(3,2, 3) The 3-line configurations giving these classes 
are of the form 



®oo 
o«o 



and therefore 

yUx(3,2,3) = x^x^AiiiX)- 

Finally, the case {Uq = U\ = U2} is related to the convoys studied in Section[8l Indeed, the formula 
follows from the density of Uq, Ui and the result that convoys are renewal processes. A more direct 
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approach follows. As there are no 3rd class particles in this case, we will use the projection into the 2 + 1 
type TASEP using only xi,X2 (or equivalently, x^ = X2). The only 2-line configuration giving classes 
(2,2,2) is 

poof 

and therefore 

^^(2, 2, 2) = Y[^xl(x2 - xi). 
Dividing by X2 - x\ and taking a limit X2 x\ gives the density Tx . □ 

8 Convoys 

The convoy phenomenon is the fact that even though each particle's speed is uniform on [-1, 1], any 
two particles have positive probability of having equal speeds. Indeed, a.s. there will be infinitely many 
particles with the same speed as any given particle. We refer to such sets of particles as convoys. Thus Z 
is partitioned in some translation invariant way into disjoint infinite convoys. 

Let Cyt = {n : Vn = Uk} denote the convoy of particle k, i.e. all particles with the same speed as k. 
We will restrict ourselves here to the study of a single convoy, though the multi-line description of the 
multi-type stationary distribution can in principle be used to understand the joint distribution of several 
convoys. 

Proof of lTheorem L8\ Partition the particles into three classes, with thresholds x = {u,u + s). The sta- 
tionary measure //x has particles of classes 1, 2, 3 with respective densities u,s,l - u - s. It is known 
that the second class particles form a renewal process. The key to the proof is (as above) to condition on 
Uo e [u,u + s] and let s ^ 0. 

Consider the two line process giving //^j and let T/,, 5yt be the counting functions of particles in the top 
and bottom lines respectively, so that T/, is the number of particles in (0, k] in the top line. We may extend 
S, r to negative k by having 5^: be minus the number of particles in {-k, 0] and similarly for Tj,. It is clear 
that {Sk}, {Tk} are random walks with {0, 1} steps with P(5<:+i - Sk = l) = u + s and P(r^.+i - Tt = I) = u. 
Let y 6 {1, 2, 3}^ denote the resulting configuration with the stationary distribution with these densities. 

The two-line collapsing procedure implies the identity 

{Vi = 2} = {5i = l,ri = 0,min5, - > 0} = {minS, - > 0}, 

k>0 k>Q 

(since 5o = = 0). Further, Vj, = 2 if and only if 5^-5 = \,Tk-Tk-i = and min;'>yt Sf-T^ = S t-Tj,. 
This suggests looking at the random walki?^: = Sk-Tk, with steps with distribution 

{u{u + e) X = \, 

U{U + S)+ UU + s X = 0, 
UU + s X = -I. 
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Having the second class particle at 1 implies that R stays positive, while its drift is 0(s). As £ ^ 
the distribution of R converges (in the product topology for sequences) to a random walk conditioned to 
stay positive for all n > with step distribution 



F(Rk^,-Rk = x) 



\uu X = +1, 

\u^ + if X = 0. 



Thus i? is a lazy simple random walk, and the only effect of u is through the probability of making a 
non-zero move. Having a second class particle at 1 does not depend on values of R„ for n <0, and this is 
also the case in the limit as e ^ 0. 

This random walk conditioned to stay positive will a.s. tend to oo as n ^ oo. Furthermore, if we take 
u = Uq then as e ^ the second class particles are exactly at k with = Uq. In particular, the convoy 
Ci is equal in law to the times of the last visits of R to any value: 

Ci= {n : m>n =^ Rm>Rn}- 

The claim that the convoys are renewal processes follows either from the corresponding fact about 
the times of last visits of R conditioned to remain positive, or from the fact that for any e > the second 
class particles form a renewal process. 

If the random walk were just a simple random walk (not lazy) then the probability of having a jump 
of length 2k + \ (as even lengths are impossible) would be p2k+i = 2"^^*+'^^(^/'). The laziness of the 
random walk implies that the distance from a particle to the next in a convoy with speed m is a sum 
of K geometric random variables with mean 1/(2mm) where P(^ = 2k + I) is as above. In particular, 
ndist = m)y<^. □ 

Example 8.1. Consider F(Uo = Ui = • • • = Un). The probability that all these speeds are in some 
infinitesimal du is 

P(t/o, ...,Unedu) = (uu)" du. 

(This can be seen easily from the corresponding density u udufor two particles and the renewal property.) 
Integrating gives 



P(C/o = ••• = [/„) = 



(2n+ 1)! 



9 Joint Distribution - ASEP 

We present two variations of our argument. The first is restricted to considering the probability that two 
adjacent particles are unswapped at large time. This event is roughly equivalent to [Uq < Ui), with some 
contribution from {Uq = Ui). 

The second variation came from an attempt to extract the complete joint distribution of two speed. 
For the ASEP it is less successful than form the TASEP, and is also conditional on a.s. existence of the 
speeds process. 
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9.1 Swap probabilities 

The key to our analysis of swap probabilities in the ASEP is to double count swaps happening until time 
t. Let R{t) be the expected number of particles j > that are swapped with at time t, i.e. 

Rit) = E#{j > : Xo{t) > Xjit)}. 

Recall the time t speed process U(t) is defined by = ML£ Define the empiric time t measure Vf 



by ^ 



The following is equivalent to the standard hydrodynamic limit theorem for the ASEP started with the 
Riemann initial condition. 

Lemma 9.1. Almost surely 2pv, converges weakly to the Lebesgue measure on R x [-p,p] 

The following simple fact is frequently useful. 
Lemma 9.2. Let Xi,X2 be topological spaces and (X{t), Y(t)), t >0 be random variables on the product 

proh (list 

space Xi X Xj- Suppose that X(t) > x and Y(t) > Y where x £ Xi and Y is an Xi-valued random 

variable. Then the joint limit also holds: iX{t), Y{t)) — (x, Y). 

The application in our case involves X(t) = Vt, which converges in probability to Lebesgue measure 
on a stripe (in the space of measures) and Y{t) = t/o(0 which tends to Uq. The conclusion implies that 
the hydrodynamic limit also holds conditioned on t/o- 

The next lemma determines the asymptotic value of R{t). 

Lemma 9.3. R(t) ~ pt/3. 

Proof. Particle has swapped with particle j > if and only if Xj(t) < Xo^t), which can be written as 

Uj(t)<Uo(t)-L 
t 

It follows that 

^ = E[v,({(x,y) : < X < Uo(t)-y})]. 

Now. lLemma 9.21 (see the subsequent discussion) shows that we can take a joint limit as Uo{t) converges 
in distribution to uniform on [-p,p] and Vt converges weakly in probability to a fixed measure which is 
l/(2p) times Lebesgue on a strip. Thus 



lim 



^Qohaix,y):0<x<U-y,y>-p}) 
2p 



XU+pf 

' 4p 3' 
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Consider now the following probability (ITheorem 1.41 shows that the two definitions are equivalent): 

Q(t) = nXoit) < Xm = P(Fo(0 < Yi(t)). (18) 
Q(t) measures the probability that particles and 1 are unswapped at time t — our present objective. 
Lemma 9.4. Q(t) is monotone decreasing in t. 

Proof. Condition on all events except those involving particles {0, 1} and denote this cr- field by !Fo,i. 
Recall that 7o,i(0 denotes the time and 1 spends next to each other up to time t and note that 7o,i(0 is 
measurable in !Fo,i. Then by (HI) we have P(Xo(0 < ^i(Ol^o,i) = ^ + pe'-'^-^'-^K Since 7o,i(0 is increasing, 
Q{t) = p + pEe'-'"-'^'^ is decreasing. □ 



Lemma 9.5. For any t we have j^Rit) = pQ(t) - p Q{t) = p + Q(t) - 1. 

Proof. Let r'l'(t) (resp. rT{t)) be the probability that at time t particle / has a larger indexed particle to its 
right (resp. left). By translation invariance these do not depend on /. R{t) is the expectation of a random 
variable which increases by one with rate p if the particle at Xo(0 + 1 has a positive index and decreases 
by one with rate p if the particle at Xo(t) - 1 has a positive index. Thus we have 

jR(t) = pr^,(t)-pr-,(t). (19) 

Consider the set A of i with a higher particle to i's right, and the set B = {n : y„(f) < Yn+i{t)}. 
By translation invariance, the density of A is r^{t), and the density of B is Q{t). There is a bijection 
between the sets, mapping / 6 A to Xi{t) 6 B. Applying the mass transport principle (see e.g. [14]), to the 
transportation of a unit mass from each / 6 A to Xi{t) 6 5 we find that r^{t) = Q{t). The same argument 
shows r~{t) = Q(t). □ 

Proof of lTheorem 1.1 2\ Combining the previous three lemmas gives that 

p/3 = limp + 2(0- 1, 

(where the limit exists due to monotonicity proved in Lemma [941) and therefore lim,^oo Q(t) = n 



9.2 Joint density 

Throughout this subsection we assume Conjecture II. 9[ Under this assumption we can talk about the 
eventual speed of a particle, and we know that for large t the empiric speed approximates the eventual 
speed. We consider the quantity 



Ra,h(t) = E 



^ l{Uj < a,Xo(t) > Xj(t)} ■l[Uo> b] 



Ra,b(t) = J]F(Uo>b, Uj < a, Xo(t) > Xj(t)) , 



Thus we ask for to have speed at least b and count particles of speed at most a that it overtakes by time 
t. This is of interest for any pair -p < a < b < p. 
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Lemma 9.6. Assume\Conjecture 1.9\holds. Then 



P , f ry-^.i A iP + a){p-b){2p + b-a) ^ 
Ra,h{t) ~ M —rr- "3^"-^ = — ^ 

J-p Jh 



4p2 Sp2 



Note: This essentially says that the contribution to i?^^ from having speed y (or in dy) and /s that 
have speed x is roughly ^t. 

Proof. Each particle moves at rate at most 1, so we have P(Xo(0 > Xj{t) < P(Poi(2f) > j). This implies 
that 

3t 

Ra,b{t) = oil) + _2 P(c/o > Z7, Uj < a, Xoit) > Xjit)) . 

The probability that any particle deviates at time t by more than s from its eventual speed is o(l). It 
follows that 

Ra,bit) = o(t) + J]F(Uo(t) > b, Ujit) < a, Xoit) > Xjit)) . 

From here on we argue as in the proof of ILemma 9.31 The hydrodynamic limit shows that Ra,h{t) is 
asymptotically close to what it would be if the speeds were independent uniform on [-p,p]: 

-Ra,b{t) = o(l) + e[i{C/o(0 > b] ■ vt{{x,y) : x e (0, Uo{t)-y\y <a,x<3}] 

= o(l) + ^ E[l{t/o(0 > b} ■ fieb {{x,y) : x e (0, Uo(t) -y),-p<y<a}] 

= o(\) + ^ E[l{Uo(t) > b] ■ (2Uo+p- a)(a+p)]. 

Simple integration completes the proof. □ 

Let Qa,h{t) be the probability of having at time t, in positions 0,1 two particles of speeds in \b, 1] and 
[-l,a] respectively: 

Qa,b{t) = ^(UyoU) > b and Uy^u) < a). 
We also let Qa,hit) be the probability of having the same speeds but exchanged: 

Qa,b(t) = P(f/Fo(o < a and t/y,(,) > b). 



Lemma 9.7. Assume Conjecture 1.9 holds. Then for any a, b, t 

d 

-TRa,h(t) = (pQa,b(t) - pQa,b{t)) 

dt 

Proof. This is an analogue of ILemma 9.51 Ra,b{t) is the expected size of the set of j's that are swapped 
with at time t (with some constraints on Uq, Uj). This set increases when has speed at least b and swaps 
with a particle of speed at most a. Using ergodicity and translation invariance, just as in ILemma 9 .51 we 
find that the expected rate at which j's are added to the set is pQa,b{t). Similarly, the expected rate at 
which elements are removed from the set is pQa,b(.t). The claim follows. □ 
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Recall that we denote by //^^ the joint distribution of Uq, U\ which we assume exists. 



Lemma 9.8. Assume Conjecture 1.9 holds. Then 

lim Qa,h{t) = m'^^KUq < -b and U\ > -a), 
lim Qa,h{t) = iP\Ui < -b and Uq > -a). 

t—>co 

Proof. Using A ^ B for A - B > 0, we have 



Qa,b{t) = nUvoit) > b and t/y,(,) < a) 

~ nUYo(t)(t) > b and t/y,«(0 < a) 
= nYo(t) < -bt and Yi(t) > 1 - at) 
= P(Xo(0 < -bt and Xi(t) > I - at) 
= P(f/o(0 < -b and Ui{t) > -a) 
^F(Uo< -b and Ui > -a) 



by convergence, 
since Xy^(f)(0 = j 
by symmetry, 
by definition 
by convergence. 



Q is dealt with similarly. 

Proof of lTheorem 1.1 3\ Combining the above lemmas and taking the limit as ? ^ oo we find that 



□ 



4p2 



dydx = pi/^\Uo < -b and Ui > -a) - pn^^\U\ < -b and Uq > -a) 



(2)^ 



where A = [-p, -b) x (-a,p]. These rectangles determine the measure pfi^^^ - pjj.^^^ in the set {(x,y) : 
-p < X < y < p} and diff'erentiating with respect to a and b gives the statement of the theorem. □ 

9.3 Equal speeds imply interaction 



Proof of \Theorem 1.14\ Since 7o,i = => Uq = Ui, it suffices to to prove P(f/o = Ui,Joj < oo) = 0. 

In the case of the TASEP the proof is very simple. From [Theorem 1.1^ we know that the probability 
that particles and 1 never swap is 1/3. On the other hand. [Theorem 1.7l implies thatP(i7o < Ui) = 1/3, 
and clearly on this event they never swap. Thus P(swap|f/o > f/i) = 1, and the result follows. 

The argument for the ASEP mirrors the above, but is more delicate with lTheorem 1.13l taking the role 
of [Theorem 1.7[ Start with 



lim 2(0 = limP(Xo(0 < X^it)) 

= limP(Xo(0 < Xiit), 7o,i < oo) + P(Xo(0 < X:(t), J = oo) 
= P(eventually Xo(t) < X:(t)) + pF{Jo,i = oo). 



(20) 
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We also have 

P(eventually Xoit) < X^it)) = P(f/o < f/i) + E [l[Uo = Ui]l[Jo,i < oo] (p + pe'^'')] . (21) 

(Compare with dH) and the discussion around it.) Combining (|20l) and (|2T]) and noting that P(7o,i = oo) = 
P(7o,i =c^,Uo = t/i) weget 

^ = P(f/o < t/,) + WiUo = f/i) + E [l[f/o = f/i]l[7o,i < cx)]pe-^«'] . (22) 
On the other hand, integrating [Theorem 1.131 gives 

= p¥{Uo < C/i) - WiUo > t/i), 

which implies 

^ = P(f/o < t/,) + ^(t/o = f/i). 

Together with (|22)) this implies 

e[uUo = f/i]l[7o,i < ^]pe-'°'] = 
and so P(t/o = Ui,Jqj < oo) = as needed. □ 

This can be extended to other particles with equal speeds. Let 7, y be the total time that particles i and 
j are in adjacent positions. 

Lemma 9.9. For any k > i, a.s. 

k = min{7 > / : Uj = Ui) => 7,-,jt = 
Consequently, in the TASEP every two particles in the same convoy swap eventually. 

Proof. Clearly this only depends on ^ - i. We proceed by induction onk - i. For k = i -\- \ this is just 
[Theorem 1.141 The key to the induction step is to show that if Uq + IJ\ then there is a transformation 
of the probability space that swaps the eventual trajectories of and 1 (and hence their speeds), keeps 
all other trajectories the same, and has finite Radon-Nikodym derivative. It follows that applying this 
transformation results in an absolutely continuous measure for the trajectories. If we assume the lemma 
for k and 1 then 

P(fc = min{j > 1 : C/^- = Uy\ and Jx,k < = 

hence by absolute continuity the result holds for k, 0. 

Recall the cr-field Tos of the trajectories of all particles except and 1. If f/o > C^i the transformation 
just eliminates all interactions between and 1 . This has the effect of exchanging their trajectories from 
some point on. Given !Fo,i, the probability of no interaction between and 1 is e"-'" ' . The Radon-Nikodym 
derivative is at most e-'" ' < oo (on 4^ Ui). 
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If Uq < Ui we define the transformation as follows: consider the first time r at which either or 
1 swaps with some other particle, and replace all interactions between and 1 by a unique interaction 
between and 1 at a time uniform on [0,t]. In the ASEP, we make this new interaction exchange 
and 1. The probability of this pattern of interactions between and 1, given To,! is pTe~^°-^, thus the 
Radon-Nikodym derivative in this case is at most e-^" ' /(pr) < oo. 

Finally, in the TASEP, since any pair of consecutive particles in a convoy a.s. swap and particles never 
unswap, it follows that all pairs eventually swap. □ 

Acknowledgment The authors wish to thank James Martin, Pablo Ferrari and Balint Virag for useful 
discussions. 
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